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Abstract 



We study in further detail a nontopological soliton model with coupling between 
quarks and mesons selected as promising in a previous study that employed the so-called 
Wigner-Seitz approximation for dense systems. Here we go beyond this approximation 
by introducing the disorder necessary to reproduce the liquid state, using the signifi- 
P5 . cant structure theory of Jhon and Eyring. We study nuclear matter, with particular 

interest in the transition to a quark plasma. The model studied is a variation of the 
chromodielectric model of Fai, Perry and Wilets, where explicit coupling to a scalar 
meson field is introduced. 
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1 Introduction 



This is the second of two papers that seek to select and develop a soliton model of nucleons 
for application to dense matter and, in particular, the transition to deconfined quark matter. 
In the previous paper, hereafter referred to as (I), we presented the motivation for our 
choice of a particular class of nontopological soliton models, the Friedberg-Lee type models 
[|l| , which have explicit quark degrees of freedom and a dynamical confinement mechanism 
resulting from a composite scalar gluon field that forms a solitonic bag in which constituent 
quarks then reside. We considered extensions of these models that include explicit meson 
degrees of freedom coupled linearly to the quarks in order to obtain a reasonable description 
of nuclear interactions. For simplicity, we have studied — and here also study — only 
models that include scalar and vector mesons, neglecting in these initial investigations any 
possible explicit effects of pions in dense matter (although the scalar meson can be considered 
an effective two-pion resonance). We avoid any double counting of hadronic degrees of 
freedom by including only nuclear constituent quarks in our calculations. Gluons enter our 
calculations only through the scalar glueball field; perturbative gluonic effects are ignored in 
nuclear matter. 

The distinguishing feature between the various models we considered in (I) is the precise 
form of the coupling between the quarks and the glueball field. We compared the various 
FL models by studying dense matter within the Wigner-Seitz approximation. In general, 
one can find parameter choices that produce reasonable results for free nucleon properties 
independent of the particular form of the quark-glueball coupling. It is in dense matter, 
where quark bags begin to touch, that the various models are distinguished. We found 
that models which have a quark-glueball coupling in accord with the dictates of the chiral 
chromodielectric model (xCD) show a behavior more in line with phenomenology. In these 
models, the quark-gluon coupling is of leading order two or greater in the glueball field, 
which is essential for the elimination of transitions to unphysical quark plasma phases at 
unrealistically low densities. Furthermore, it was found that coupling the quarks to a scalar 
meson field ensures saturation. The quark-meson coupling is taken to be independent of the 
glueball field within the mean field approximation to avoid unphysical transitions in dense 
matter, as was detailed in (I). 

Here we wish to further study the model selected in (I) by going beyond the relatively 
rough approximations used there in modeling the liquid state. The Wigner-Seitz approxi- 
mation consists in assuming that each nucleon is confined by interactions with its nearest 
neighbors to a given volume, equal to the inverse of the baryon density, known as the Wigner- 
Seitz cell. For a solid cubic lattice, for example, the Wigner-Seitz cell is a cube. Here, we 
are interested not in the solid but rather the liquid state, and so the usual choice is to take 
the Wigner-Seitz cell to be a sphere in the hope that one thereby better models the dis- 
order of a liquid. This, however, is clearly not enough: if we look at models of the liquid 
state used by physical chemists in order to describe molecular liquids (models that passed 
out of use several decades ago after the development of large-scale computers enabling the 
use of molecular dynamics and Monte Carlo techniques), the Wigner-Seitz approximation 



employed in (I) corresponds to the cell model of Lennard- Jones and Devonshire 0, which 
does much better at reproducing the solid state than the liquid state ^. Instead, we shall 
employ a refinement of the cell model — namely, the Significant Structure Theory of Jhon 
and Eyring — which introduces holes into the system in order to account for the disorder 
present in liquids. 

As a matter of fact, in (I) we assumed that the Wigner-Seitz cell is simply a sort of 
"average snapshot" of the nuclear medium felt by the quarks inside an otherwise freely 
moving nucleon. Thus we proceeded by subtracting away energy due to spurious center of 
mass motion (due to the fact that the nucleon was not constructed by putting the quarks in 
a good momentum state), and then took the kinetic energy of the system to be that of a free 
Fermi gas. Clearly, this approximation can only be justified at low densities. As the density 
increases the motion of an individual nucleon is affected by the medium, and this leads us 
to consider the Wigner-Seitz cell not just as a boundary upon the quark wave functions that 
build up a nucleon, but also as a restriction upon the motion of the nucleon itself. This leads 
to the considerations of the previous paragraph, which shall be further developed in Sec. |^. 

The outline of the paper is as follows. The nontopological soltion model used is reviewed 
in Sec. ^. In Sec. ^ we discuss various attempts to model a soliton liquid, then motivate 
and introduce the particular model based on siginificant liquid structures that we shall use 
here. In Sec. § we present the resulting equations of state for nuclear matter and discuss the 
transition to quark matter. A general summary and discussion of the two papers is given in 
Sec. g 



2 The Model 

The nontopological soliton model we study here is based upon the chiral chromodielectric 
model of Fai, Perry and Wilets 0. In its full version, the model contains quark and gluon 
degrees of freedom. A scalar glueball field a couples to the quarks, and colored gluons A^are 
treated perturbatively. The scalar field provides absolute confinement of both quarks and 
gluons and gives consituents a mass. Meson exchange is surely present in this model, but 
for simplicity we alter the original xCD by dropping the gluon field A" and ignoring sea 
quarks. Instead, as in quark-meson coupling models, we introduce a scalar meson (p. The 
vector meson V^, which provides repulsion in quantum hadrodynamics, is not necessary here 
since the soliton structure provides repulsion between nucleons, and so for simplicity we set 
V^ = 0. We assume the scalar meson couples linearly to the quarks and take the quark- meson 
vertex to be independent of a. The Lagrangian density for our model is 

C = 2p[zYd,-g{a)-gs(P]iJ + ^d^ad''a-U{a) 

+]^d,<pd^<P-\ml<p'-\F,,F^\ (1) 



where 



U{<y) = 1^^ + |a^ + ^_^' + B . (2) 



We have set the current quark masses to zero. The parameters are as chosen in (I). The 
constants of the potential are a = 50fm~^, b = — 1300fm^^ and c = 10^, so that U{cr) has a 
local minimum at o" = and a global minimum at a = a^ = 0.285fm~^, the vacuum value. 
The mass of the glueball excitation associated with the a field is mcB = \/U"{av) = 1.82GeV 
and the bag constant is -B = 46.6MeV/fm^. The quark-cr coupling is 



g{a) = g^a,,[ 
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where we choose the chromodielectric function nicx) to be 
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In the following we take g^ = 3 and k^ = .1. 



We solve the Euler-Lagrange equations in the mean field approximation, replacing the 
glueball field a by the classical soliton solution a{r) and the meson field by its expectation 
value in the nuclear medium < >= 0o- The resulting equations for the quark and the 
scalar soliton field are solved in a Wigner-Seitz cell of radius R by implementing boundary 
conditions based upon Bloch's theorem, as detailed in (I). The quark spinor in the lowest 
band is assumed to be an s-state 
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(5) 



and we make the simplifying assumption of identifying the bottom of the lowest band by the 
demand that the derivative of the upper component of the Dirac function disappears at i?, 
and the top of that band by the demand that the value of the upper component is zero at 
R [|]. The resulting equations for the spinor components are 



duk 
dr 



+ [5'(o") - 9s<pQ + Cfc] Vk = ^ 



dvk '^Vk 
dr r 



[fi'(o-) - 9s4>0 - Cfc] Mfc = . 

The equation for the soliton field is 

-V'a + U'{a)+g'{a)p,{r) = 0. 
The quark density Pg and the quark scalar density ps are given by 
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47rA;3/3 Jq 

where the band is filled up to k. The quark functions are normalized to unity in the Wigner- 
Seitz cell. The boundary conditions for the sohton field are cr'(O) = cr'(i?) = 0. The boundary 



conditions for the quark functions at the origin are given by u{0) = Uq and f (0) = 0, where 
Uq is determined by the normahzation condition 

R 

Anr'^driuirf + v{rf) = 1 . (11) 



The boundary conditions ai r = R are given by u'^{R) = and Vf,{R) = for the bottom 
of the lowest band, and Ut{R) = for the top of this band. Using these equations we 
can solve for the corresponding eb and e*. We assume the tight-binding dispersion relation 
ek = el + (et — e^,) sin^(7rs/2), with s = k/kt, and that the band is filled right to the top — 
that is, k = kt. With this dispersion relation and filling, the nucleon energy is given by 



f-i _ r _ /ttcx > ^h _ rl 

-a\rf + U{a) 



En = ^riq J ds s^ Itb + {tt - et) sin^ Ty j | + ^ Anr'^dr 



(12) 



In order to correct for the spurious center of mass motion in the Wigner-Seitz cell the nucleon 
mass at rest is taken to be 



Mn = ^JEl- < P|„ >ws , (13) 

where < P"^^ >ws= ng < p^ >ws is the sum of the expectation values of the squares of the 
momenta of the ng=3 quarks. 

At low density the band width vanishes and the quarks are confined in separate bags. 
Then we can assume the individual nucleons move around as a gas of fermions with effective 
mass Mat given by Eq. (0), so the nucleon energy is Ej^ = J Mj^ + k^. The total energy 
density at nuclear density pb is thus 



7 



kp 



^9 = J^l' dk^Ml + P + -mlcg, (14) 

where 7 = 4 is the spin-isospin degeneracy of the nucleons. The Fermi momentum of the 
nucleons is related to the baryon density through the relation 

The total energy per baryon is given by Eg = Sg/ps- We have used the label g to indicate 
that these expressions correspond to a gas-like phase. The constant scalar meson field 0o is 
determined by the thermodynamic demand of minimizing S: 

(900 

Our mean field equations are similar to those of quantum hadrodynamics, the difference here 
being that the nucleon now has structure and thus the meson field couples to the nucleon 
through its quarks. 

Let us review the treatment of dense matter in (I). We started with the so-called Wigner- 
Seitz approximation, which is often used in soliton calculations, since it is the simplest 
picture of dense matter available. In this approximation, each soliton is confined to a unit 



cell, which we can view from two different perspectives. If we choose periodic conditions at 
the cell boundary, we have the usual Bloch approach to the solid state: the quarks play the 
role of the electrons and the scalar glueball field takes the role of the ions in the usual crystal 
formulation. If instead we want to model the liquid state, we need to somehow introduce some 
disorder into the system. From this view, then, we take the Wigner-Seitz approximation as 
a sort of averaging over the rest of the system: for the purposes of constructing an individual 
nucleon, we ignore the motion of the nucleons, instead adding "by hand" the kinetic energy of 
a free gas to the system. This is justified to the extent that each nucleon's motion describes 
slow degrees of freedom, whereas the constituents are fast degrees of freedom that react 
essentially instantaneously to changes in the relative arrangement of the nucleons — that is, 
if the Born-Oppenheimer approximation is valid. 

Thus in (I) we calculated the equation of state of nuclear matter within the second 
scheme. The mass of the nucleon was calculated as a function of the radius of the (spherical) 
Wigner-Seitz cell, for which it was then necessary to subtract away spurious center-of-mass 
motion. This was done using approximate relations discussed in detail in [^]. Note that 
within the first scheme, the quarks are not assumed to be in a state of good momentum, and 
this kinetic energy is not spurious. 



3 Significant Structure Theory 

In (I) we used the Wigner-Seitz approximation to determine the effective nucleon mass 
in nuclear matter, but then gave the nucleons the kinetic energy of a Fermi gas. In this 
approximation, the quarks feel the nuclear medium, but the nucleons themselves do not. 
That is, we assume the quarks adjust instantly to the medium, forming 3g collective states 
that move relatively slowly and essentially freely through the medium. The interactions 
between nucleons occur only indirectly through the effective mass and the mesonic mean 
field. Clearly, such an approximation cannot be accurate at high densities, when the finite 
size of the nucleon becomes important, for nucleons will not then move freely. Instead, we 
need to model the liquid state, where any individual nucleon will range about the system over 
long time scales, but will be localized on shorter time scales. One would still like to approach 
the problem using the Wigner-Seitz approximation, insisting now that the nucleon also feels 
the medium and does not leave the cell. Clearly, this is a very restrictive assumption, and 
physical chemists long ago understood that this corresponds more closely to the solid than 
the liquid state. Nevertheless, this will provide a starting point for the model of the liquid 
state we shall adopt in the following, so let us pursue this approach further. 

Ideally, we would like to allow our nucleon to rattle about in its Wigner-Seitz cell in order 
to extract a potential. This would entail dropping the assumptions of spherical symmetry 
for the bag and the use of only s-wave quark states. Instead, we shall attempt to model the 
nucleon's motion at high density as follows. First, we assume that the motion is harmonic 
— that is, the center of mass Rrm of the nucleon is never far from the center of the WS cell 



Ro = 0. Now, the average of a harmonic potential in the ground state is 

{V) = ^Im^^^RL) = \uJN, (17) 

from which we find the natural frequency uo^ as a function of the effective nucleon mass 
Mat and the mean square of the center of mass coordinate. Next, we identify Rem as the 
center of the quark distribution in the cell, ignoring thereby any motion of the soliton bag 
a. (Viewing Fig. (3) of (I), we see that the soliton bag begins to be "squeezed" by the WS 
boundary at i? ~ 1.5fm, so that at least for densities higher than this we might argue that 
the bag is essentially fixed.) Taking Rem to be the center of mass coordinate of three quarks, 
we have 
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Now we identify the last average with the mean square charge radius of the nucleon 
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This is clearly only a rough estimate of the center of mass motion of the nucleon. The 
nucleon energy in the solid is now 

e{^) = M^ + L^, with o;^ = --l-_. (20) 

2 2Mjv(i?^m) 

This approximation corresponds to subtracting away spurious kinetic energy from the soliton 
energy Eq. (p^, namely, 

E,, = E^-E^j^^ = ^Ml + {P^Jws -MM-\io (21) 

\Pcm)wS - 
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Approximating (-Rem) by \{'r'^)ws can only be valid at high density. At low density, this 
surely breaks down for then the quarks cannot reach the WS boundary unless the bag 
itself is allowed to move. Moreover, the present approximation corresponds to treating the 
soliton matter as an Einstein solid (we have implicitly averaged over the Bloch momenta K 
corresponding to the lattice of nucleons: a more accurate treatment of the solid would put 
the three quarks in a state of good K = ki + k2 + ka and find a frequency a;7v(K) that is a 
function of the Bloch momentum). The total energy density in this "solid" phase is 

£, = pb (^Mn + -c^Tv) + -ml(Pl. (22) 

The constant scalar meson field (po is again determined by the thermodynamic demand of 
minimizing S. Both M^ and u depend implicitly upon 0o, so that this equation must be 
solved iteratively in conjunction with those for a and ip. 



This gives the equation of state for sohd nuclear matter, which is not of much interest in 
itself. However, this is useful for building a model of the hquid state at high density based 
upon significant structure theory |^. The essential idea is to isolate those configuations 
that make the significant contribution to the partition function. Based upon experimental 
observations of molecular liquids, the liquid state is viewed as a close-packed lattice with holes 
present that destroy any long-range order. The volume of the system increases by increasing 
the number of holes, with the volume of each hole equal to the average volume occupied by 
a close-packed molecule, since this balances the competing demands for greater entropy and 
lower energy. A molecule neighboring a hole can move into the vacancy, creating a new hole 
at the site it left. Each hole thus replaces three vibrational with three translational degrees 
of freedom. Thus the liquid is represented as a combination of molecules with solid-like 
properties (those next to filled sites) and gas-like properties (those next to vacant sites). 

Now consider this model of the liquid state applied to dense solitonic matter. (Such an 
application has been studied previously for dense skyrmion matter in Ref. [0.) First, we 
note that our assumptions of a spherical Wigner-Seitz cell, which reflects a sort of average 
over nearest neighbor positions, and the Einstein approximation that ujn is independent of 
the Bloch momentum K, which ignores long-range correlated vibrations of the nucleons — 
approximations that would be rather severe if we truly wished to model a solid — are instead 
appropriate for the present application. We need only add the holes to ensure the disorder 
corresponding to a liquid state. So let Vi be the total volume of the liquid and v be the 
volume of each cell (occupied or unoccupied). If there are A^ nucleons and Nh holes, then 
Vs = Nv is the total volume of the occupied cells and Vg = Vi — Vg = NhV is the total volume 
of the holes. On average, a nucleon will encounter a neighbor on a fraction j^^j^ = ^ of its 

trips and a hole on a fraction 1 — ^ = y- of its trips. Thus there are 3iVy- solid-like and 

3iV-rf gas-like degrees of freedom, and the liquid partition function is 



Vi 



,^T7 ^^(1-w) 



Zi = Zs "' Zg "' . (23) 

This results in the following nucleon energy in the liquid state: 

Ei>-^E^\-) + ^mn,). (24) 

where n^ = -rr- = ^ is the "density" of the gas-like part of the system, which from the last 

equality is equal to the baryon density ps- The average solid-like nucleon energy E^ {v) 
is given by Eq. (pO[) , with Mj^ and ojn depending upon v through the Wigner-Seitz radius 
R = (3f/47r)^/'^. The average gas-like nucleon energy is instead 

E^^\n,) = ^ fj d'k^MUvi) + P, (25) 

with kg = [Qn'^ng/'^Y/^ = (Gtc'^Pb/iY^^ and vi = Vi/N = l/ug. The total energy density of 
the system in the liquid phase is then 

Si{v) = pIvE^^\v) + pb{1 - Pbv0^\pb) + \rnl<Pl - \rnlVl (26) 



Once again, (J)q is determined by minimizing Si, and so for nonzero Qs we must solve anew the 
equations for a and ip consistent with the mean field value 0o(^)Pb) determined from this 
new liquid equation of state. In Eq. (pUf), the cell volume v is to be taken as a parameter. 
Note that the WS cell volume is no longer 1/pb when holes are present. Instead, we can 
define a new radius Ri = {S/AttpbY^^ which is half the average spacing between nucleons in 
the liquid state. 



There is some question as to whether one should take the effective nucleon mass in ( p5D 
to be Mjv(t>) or Mj^{vi) = M^iPs^)- Since a single nucleon can have both gas-like and solid- 
like degrees of freedom, it might be more consistent to use the former choice in calculating 
the gas-like part of the nucleon energy. However, there seems little reason to insist that the 
inertial masses corresponding to the solid-like harmonic motion and the gas-like translational 
motion are the same. Moreover, the latter choice guarantees that our liquid EOS reduces to 
the Fermi gas EOS given by ([T^ ) at low density. Thus our equation of state (p6|) interpolates 



smoothly between high-density and low-density pictures of the liquid state. 



4 Results 

4.1 Nuclear matter 

In Figs. 1-3 we show the energy per baryon of xCD solitonic nuclear matter for three different 
values of the quark-meson coupling. In each graph we display the "gas" equation of state 
(|14D and a set of significant liquid structure curves characterized by various choices of the 
WS cell volume v = 4nR^ in (|26|). In addition, the upper curve labelled "solid" in these 
figures is given by ([T^) with the condition Ug = 1/v — or, equivalently, Ri = R — which 
implies that the number of holes is zero, so that the nuclear motion is entirely solid-like. 
Note that our approximations in estimating uj^ break down at low density. In particular, 
the "solid" curve is not to be trusted above R ~ 1.5fm: here, the bag is no longer squeezed 
by the WS boundary and its motion probably can no longer be ignored. 

As f — ^ 0, the significant structure curve approaches the "gas" curve, which certainly 
underestimates the energy per nucleon. Without empirical evidence of a transition from 
liquid to sohd nuclear matter, we are unable to fix v directly. (Some models predict such a 
transition; however, it is more widely believed that the transition to the quark-gluon plasma 
will preclude a transition to solid nuclear matter.) Here, we can just treat f as a parameter 
in choosing the EOS that best fits the empirical data at the saturation point. We find that, 
for the parameter sets we studied, the curves that best fit the empirical EOS for nuclear 
matter are given by the values gg = 2,R = 0.4fm and Qs = 3,R = 0.7fm. The values of 
the Fermi momentum, binding energy, and compression modulus at the saturation point 
are fc, = 1.06fm-^ E^ = 22MeV and K = 1082MeV for the former, and fc, = 0.95fm-\ 
Eg = 24MeV and K = 6 71 Me V for the the latter curve. These are to be compared with 
the empirical values kg = 1.36fm~^, Eg = 16MeV and K ^ 200MeV. In particular, we 
see a significant improvement in the value of the compression modulus with respect to the 
calculation of (I), although it is still rather high with respect to empirical values. 

9 



4.2 Quark matter 

In the high density hmit, the preferred phase in the nontopological sohton models we are 
studying is a uniform plasma characterized by the solution a = 0. That is, the soliton bags 
disappear, and one is left with a quark gas. The energy density of the quark plasma is 

where the Fermi momentum is related to the baryon density pB as kp = (Qt^'^Pb/iY — 
the degeneracy of the quark gas is Ugj and the quark density is nq=3 times the baryon 
density. The bag constant is -B = A6.6MeV/ fm^ for our choice of parameters. In contrast 
to our assumption about the solitonic phase, the energy density of the quark gas is altered 
significantly by perturbative gluonic contributions, especially at higher densities. In this 
phase, the xCD model in the mean field approximation is equivalent to perturbative QCD. 



As shown in P2[, for example, adding the lowest order gluonic corrections then gives 



3k% {^ la. a? 






6.79 + 2 In f^ 



B, (2^ 



where the leading-logarithm expression for the strong coupling constant is 

We take the QCD scale parameter to be A ~ 180-200MeV. 

In Fig. 4 we show the two nuclear matter curves selected above along with the quark- 
gluon plasma EOS given here. In addition, we show an "empirical" nuclear matter EOS 
given by 

^B-^(|-l)%Mr-E„ (30) 

with mJ^"'^ - Es = lieOMeV and K = 200MeV. (The bag constant B appearing in the 
quark-gluon plasma energy has been set in fitting the free soliton mass and rms radius, and 
thus for purposes of comparison we take the low density limit of the "empirical" curve to 
coincide with the free soliton mass.) Note that the Qs = 2 shown here has a transition to 
the solid state at kp = (97r)^/^/2i? = 3.8fm~^ and the f^^ = 3 at 2.2fm~^, both well past 
the transition points to the quark-gluon plasma. Clearly, the saturation points of the model 
curves occur at densities that are too low. Nevertheless, the very fact that we find qualitative 
agreement with the empirical nuclear matter EOS is quite encouraging. 



5 Conclusions and Outlook 

In this paper we have developed an improved modeling of the liquid state of solitonic matter 
based upon the significant structure model used in physical chemistry. We have applied this 

10 



to the study of nuclear matter within a nontopological sohton model with explicit quark 
degrees of freedom. In particular, initial studies in (I) indicated that among several such 
similar models, the chiral chromodielectric model gave results most in line with empirical 
expectations. When a scalar meson is allowed to couple to the quarks, saturation can 
be achieved. The calculations in (I) were based upon a modeling of nuclear matter that 
assumes the solitons move essentially freely through the medium. This clearly will break 
down at and above nuclear saturation density, where the size of the individual solitons and 
the spacing between solitons in the medium become comparable. The significant structure 
model of the liquid state used in the present paper, however, is designed for densities near 
the transition to the solid state. We have thus developed here an equation of state that 
interpolates between models designed for low and high densities. It is encouraging to find 
that without changing the free nucleon properties we can adjust the quark-meson coupling 
to reproduce the saturation point of nuclear matter found in (I). Indeed, we find that the 
compression modulus is improved significantly when calculated using the significant liquid 
structure model. 

With the xCD model we are able to treat the transition to the quark-gluon plasma 
consistently. This is the point of developing this model: confinement is dynamical. Thus 
the parameters governing the nuclear and plasma phases are in principle the same. Clearly, 
the results we have found here, while qualitatively in agreement with empirical estimates, 
are not quantitatively useful. In particular, the saturation point of nuclear matter occurs at 
far too low a density for the parameters used in our calculation. Moreover, the mass of the 
free nucleon is too high. We have not tried too hard to improve our results by adjusting the 
parameters of the model. Probably one can find a better choice of parameters than those 
for which we have done the calculations here, but we must emphasize that there seems to 
be a limit in how well one can do. This can be seen already in studying the free nucleon. 
There, we are unable to find a parameter set that reproduces the nucleon mass and rms 
radius exactly while still giving reasonable values for the glueball and bag constant. We 
compromised here by fitting the nucleon rms radius well while keeping the glueball mass and 
bag constant in their accepted ranges. This resulted in a free nucleon mass that was too 
high. 

In trying to reproduce the empirical saturation point of nuclear matter, however, we 
found a general scaling phenomenon. Keeping the saturation energy roughly correct, it 
seems difficult to change the parameters in such a way as to get the correct density expect 
by a rough overall scaling of all predicted quantities. Thus getting the correct saturation 
density entails lowering the rms radius along with a corresponding increase in the mass of the 
nucleon. This is suggestive. As a matter of fact, the quark-meson model we developed here 
is better suited for nuclear matter than for isolated nucleons. This is because an isolated 
nucleon will surround itself with a pion cloud, whereas (unless pion condensation occurs) 
the effects of pions in nuclear matter are likely to be small if we already have scalar mesons. 
Thus we can consider the model used here as actually better used for describing the quark 
core of the nucleon. The addition of the pion will allow the quark core to shrink and lower 
the energy of the free nucleon with respect to that of the quarks. With the shrinking of 
the quark core, one can expect a corresponding decrease in the volume per nucleon at the 
saturation point. 

11 



Thus we view the addition of pions as an essential improvement to our modeL This, of 
course, was always obvious: we certainly cannot hope to reproduce the long-range part of the 
nucleon-nucleon interaction without pions. Our results indicate that the presence of pions 
is necessary in order to reproduce the structure of the nucleon as well. This comes as little 
surprise. There are of course other improvements that can be made upon our calculations, 
such as a better handling of the Bloch boundary conditions and quark wave functions 0] and 
an improved treatment of the corrections due to spurious center of mass motion. Perturbative 
corrections due to gluons and mesons about the mean field should eventually be considered 
as well. 

Having said (that is, written) this about improving the model, we should not lose sight of 
the original goal of our work. What we wanted to do first of all was see if we could distinguish 
among the various nontopological soliton models on the market by studying dense matter. 
To this we can answer that the chiral chromodielectric model appears to be more in line 
with empirical expectations. Then, we wanted to develop as simple a model as possible that 
could give a rough qualitative fit to both free nucleon and dense nuclear matter properties, 
thus providing a reasonable starting point for more sophisticated models that can provide 
truly quantitative predictions. To this we can also answer that the chiral chromodielectric 
model, when modified according to the local uniform approximation by the addition of a 
scalar meson field, would appear to be such a model. In fact, this model not only gives a 
good qualitative and rough quantitative fit to the empirical equation of state, it also predicts 
an increase of the nuclear rms radius in the nuclear medium, a result that is in accord with 
the EMC effect. Moreover, the results presented here seem to indicate that with the addition 
of pions we would have good possibilities of obtaining a quantitatively accurate fit to single 
nucleon properties and to the empirical nuclear matter equation of state. This would provide 
a reliable estimate of the bag constant and therefore a consistent and accurate treatment of 
the transition from nuclear to quark-gluon matter as a true transition between two phases 
of a single model. 
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Figure 1: Energy per baryon as a function of average liquid radius Ri = {S/AttpbY^^ for 
various values of the Wigner-Seitz radius R = (3v/47r)^/^ for quark-meson coupling gg = 1. 
Also shown are the curves for the "gas" limit R ^ and the "solid" limit Ri = R. 
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Figure 2: As in Fig. 1, for the choice of quark-meson couphng Qg = 2. 
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Figure 3: As in Figs. 1 and 2, but with g^ = 3. 
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Figure 4: Quark-gluon plasma (QGP) and nuclear matter energy per baryon number as 
function of the Fermi momentum. The curves are as described in the text. 



2000.0 



1800.0 



1600.0 






W 



1400.0 



1200.0 



1000.0 



g=2R=0.4 
g=3R=0.7 




QGP A=200 
QGPA=180 



"empirical eos' 



0.5 1.0 1.5 2.0 2.5 

Kp[fm] 



3.0 



3.5 



4.0 



17 



